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We show that dipolar bosons and fermions confined in a quasi-one-dimensional ring trap exhibit a 
rich variety of states because their interaction is inhomogeneous. For purely repulsive interactions, 
with increasing strength of the dipolar coupling there is a crossover from a gas-like state to an 
inhomogeneous crystal-like one. For small enough angles between the dipoles and the plane of the 
ring, there are regions with attractive interactions, and clustered states can form. 
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Atoms or molecules with electric or magnetic dipole 
moments offer unique opportunities for exploring 
strongly interacting few- and many-body quantum sys- 
tems. Much progress has been made recently in realiz- 
ing such systems experimentally The anisotropy of 
the dipole-dipole interaction is an important resource in 
using cold atoms or polar molecules to simulate elusive 
quantum states in condensed matter physics (see for ex- 
ample the review articles Dipolar quantum gases 
in lower dimensions have received considerable attention, 
mainly because the collisional instability toward "head- 
to-tail" alignment of the dipoles is strongly suppressed 
From a theoretical perspective, interesting aspects 
are the possibility of creating, e.g.. a p-wave Fermi su- 
perfluid in two dimensions [4| or phases with Luttinger- 
liquid-like behavior in one dimension (ID) d, 0|. 

In this article, we consider a feature of the dipole-dipole 
interaction that has so far received little attention: In 
curved lower-dimensional geometries, the two-body in- 
teraction becomes inhomogeneous i.e., it depends not 
only on the relative coordinate, but also on the center of 
mass of the two dipoles. This provides a way of creat- 
ing inhomogeneous two-body interactions that valuably 
supplement other methods, such as the optical control of 
magnetic Feshbach resonances 

The case we study here is particles confined in a quasi- 
ID ring a system which can be realized experimen- 
tally with ultracold atoms [l^, 11 1 and which is also 
of interest in connection with semiconductor nanostruc- 
tures [l2| . In this model system, sketched in Fig. [1] the 
degree of inhomogeneity can be tuned via the orientation 
of the dipoles. Here, we calculate ground-state properties 
of a few dipolar bosons and fermions by applying the nu- 
merically exact multi-configuration Hartree method [l^ . 
For dipolar tilt angles such that the interaction is repul- 
sive, we identify weakly and strongly interacting gas-like 
states. Moreover, for large enough tilt angles, there are 
regions with attractive interactions, and clustered states 
may be formed. 
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Figure 1: (color online) Sketch of the ring-shaped trap in 
the xy plane. The dipole moments d = d (sin a, 0, cos a) are 
aligned by an external field E in the ajz-plane. 



Hamiltonian — Let us consider a system of N iden- 
tical particles (bosons or fermions) of mass m confined in 
a ring-shaped trap of radius R in the xy-plane (Fig. [1]). 
The particles have a dipole moment d = d (sin a, 0, cos a) 
aligned in the xz-plane by an external field, at an angle 
a to the z-axis. The interaction between two dipoles 
is V{r) — D^{1 — 3 cos^ 9 rd) /r^ , where r is the separa- 
tion of the dipoles, 9rd is the angle between r and d, 
and = (P/Aireo for electric dipoles and cPfio/Air for 
magnetic ones. It is convenient to introduce the length 
rd = 2mD^/h^ as a measure of the dipolar interaction. 
We focus on the limit of a tightly confining ring po- 
tential, for which the transverse motion is frozen in the 
lowest-energy mode, whose spatial extent a± =y/h/mijj± 
is much smaller than R Integrating out the trans- 
verse degrees of freedom, one arrives at an effective ID 
Hamiltonian 



H = - 



■ V — 



2mR^ ^ del 



i<j 



(1) 



where the angle 9i specifies the position of particle i on 
the ring. For R ^ a±, the effective interaction takes 
the form ViD(6'i,6'2) = Vcm (6) K-ei(i?), where 0^(91 + 
6*2)72 is the center-of-mass (CM) angle of the two dipoles 
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Figure 2: (color online) Left: Coupling strength Vcm(6) for 
a — (dashed), O.lQn ~ Qc (dotted) and n/2 (solid). Right: 
Relative-coordinate potential Vrci(i?) for R/a±_ — 10, includ- 
ing the long-range asymptotics (dashed). 

and I? = 01 — 6*2 is the relative angle. In terms of the 
variable s ~ 2i?| sin(i^/2)]/a_L, the dependence on the 
relative angle is given by [7| 

^■ci(i9) = \/2^(l + s2)e"'/2crfc(s/V2) - 2s, (2) 

and VbM(e) = D2(i_3sin2 asin^ e)/4a\. Thus the CM 
potential (shown in Fig. [J) has minima at = ±7r/2, 
which become more pronounced as a increases. For 
a > ttc = arcsin(l/-\/3) ~ 0.1967r, the potential acquires 
attractive regions. 

Homogeneous interactions — We first consider the 
case of dipoles aligned perpendicular to the plane of 
the ring (a = 0), for which the interaction is positive 
and isotropic, and draw some qualitative conclusions be- 
fore proceeding to the inhomogeneous case. For — )- 0, 
Vreiii}) remains finite because the transverse average cuts 
off the divergence of the 3D interaction. Thus for 
weak enough dipolar repulsion (specified below), the in- 
teraction behaves like a potential with range ^ a± [l^. 
For bosonic dipoles, the problem then approximately 
maps to the Lieb-Liniger model [l^ on a ring with a 
contact interaction gS{xi — X2), where Xi = R6i. For 

^ aj_, 5 is the zero- momentum Fourier transform of 
the potential, and g = D^/a]_. In the thermodynamic 
limit iV, i? — >■ 00 (keeping the density n ~ N /2'kR fixed), 
the properties of the Lieb-Liniger model are determined 
by the dimensionless parameter 7 = gm/h^n. For 7^1 
the system is a weakly interacting Bose gas, while for 
7 ^ 1 it may be mapped to a noninteracting Fermi 
gas [13. 

For separations r 3> ax the effective ID interaction 
behaves as l/r"^. For sufficiently strong dipolar interac- 
tion, this tail becomes important and the system goes 
over into a crystal-like state, where the dipoles local- 
ize as an equally spaced lattice, with energy Eq ~ 
{ND^/16R^)J2^-^^ |sin[z/7r/iV]r^ The amplitude A of 
the zero-point motion may be estimated to be of order 
(h/mujy/^, where u = {12D^n^ /my/^ is the oscillation 
frequency of a single atom, all others being held fixed. 
Dipoles are well localized if A ^ i.e., when the 

dipole length is large compared with the inter-particle 
spacing, {nra)^^^ > 1. 
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Figure 3: (color online) Homogeneous system {a — 0): 
Ground-state energy E{nrd) for A'' = 4 dipoles. For compar- 
ison, the plot includes the Lieb-Liniger results for the Bose 
gas, Hartree-Fock for fermions, and the classical limit Ec- 



Let us now quantify the above considerations. In 
Fig. [31 we plot the energy of TV = 4 bosonic and 
fermionic dipoles as a function of interaction strength 
for R/a± = 10. We see that for nr^, ^ 1, the bosons 
are indeed well described by the Lieb-Liniger model with 
g = D^/a\. In particular, for 7 = rd/lna?^ ^ 1, 
the bosons form a weakly interacting gas, whereas for 
1 < 7 <C (naj^)~^ the bosons tend to fermionize and the 
energy increases less rapidly with r^. For the parameters 
adopted here, 7 = 1 at nvd ~ 2{na±^'^ w 0.01. Note 
that, since (naj^)^ = nrd/2^, the fermionized regime is 
fully separated from the crystal-like one only for strong 
confinement, na±^ <C 1. For fermions, interaction effects 
are especially small for nr^ < 1 because exchange contri- 
butions cancel the direct ones. We also plot the Hartree- 
Fock energy for comparison. 

Inhomogeneous interactions — For dipoles aligned 
perpendicular to the plane of the ring, as we discussed 
above, the system is rotationally invariant about the nor- 
mal to the plane, and the particle density is uniform. 
However, for a 7^ the interaction is smallest for CM 
angles <d = ±7r/2 and particles tend to lie closer to 
9 = ±7r/2 than to = or tt. For nrd ^ 1, the kinetic 
energy becomes unimportant, and the configuration is 
well approximated by particles localized at angles 9i that 
minimize the interaction energy. 

The upper row in Fig. 2] shows the particle density 
p{6) around the ring for TV = 4 bosons (left) and = 3 
fermions (right) for a dipolar tilt angle a ~ O.IOtt slightly 
below the threshold value ac- For small couplings nrd, 
the density (green line) reflects the Bose-gas-like state, 
that gradually develops an inhomogeneity when the in- 
teraction is increased to higher values of nrd (blue lines) : 
The initially homogeneous bosonic ground state devel- 
ops four distinct peaks that demonstrate the localization 
of the particles near their classical equilibrium positions. 
The same trend is observed for fermions, as is shown in 
the right panels of Fig. S] for = 3. For that case, the 
ground state is a superposition of two states, one of which 
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Figure 4: (color online) Density profiles p{9) for = 4 bosons 
(left) and A'' = 3 fermions (right) for a — O.IQtt {upper- 
most row) and a = 0.227r {middle row) for different cou- 
pling strengths. Bottom: Schematic overview of the different 
ground states for bosons (for fermions, see text). The dashed 
lines indicate the boundaries between these approximate re- 
gions. On the curved dotted line, the ground-state energy is 
zero. 



has two particles near 9 = tt/2 and one near — 7r/2, and 
the other in which the populations are reversed. As a 
result, the density profile has 2N = 6 peaks rather than 
N. 

For a > ttc, the interaction is attractive for Q = 
±{n/2 ± AQ) and, for small Aa = a ~ ac, AQ ~ 
2^^^\/ Aa. Then clustered states can be formed, as is 
illustrated in the middle row of Fig. [H With increas- 
ing interaction strength, the density distribution shows a 
transition from the homogeneous state to one with two 
distinct maxima on opposite sides of the ring. This in- 
dicates a clustering of the particles due to the attractive 
part of the interaction, as we further discuss below. The 
lowest diagram of Fig. [4] summarizes schematically the 
different ground states for bosons as a function of cou- 
pling strength and tilt angle, with the dashed lines indi- 
cating the approximate boundaries between the different 
regimes. The vertical lines corresponds to the values of 
nvd for which the energy of the Tonks gas equals that 
of the weakly interacting Bose gas or that of a classical 
crystal. For fermions, the diagram's structure is similar, 
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Figure 5: (color online) Ground-state energy E{nrd) for 

N = Z fermions ( — , a = 0.227r); iV = 4 bosons ( , 

a = 0.27r, 0.227r, and 7r/2 from top to bottom). 



but with the gas-like regions replaced by a single one — a 
weakly interacting Fermi gas — and clustering occurs at 
higher nr^. 

To investigate the clustered states further, we plot in 
Fig. [5] the ground-state energy of bosons and fermions as 
a function of the coupling strength for values of a greater 
than Uc- The coupling required to produce a bound state 
is greater for fermions than for bosons for two reasons: 
first, the Pauli principle leads to a greater kinetic energy 
for fermions and, second, the exchange hole for fermions 
reduces the magnitude of the interaction energy. 

With increasing dipole moment, the energy decreases 
for all a > at sufficiently high nr^ for both fermions 
and bosons. For bosons and for a < the en- 

ergy increases monotonically as a function of nr^, for 
Uc < a < o.c2 = arcsin \f2ji w 0.3047r it exhibits a max- 
imum, while for a > ac2 it decreases monotonically. The 
reason for the change in behavior at a = ac2 is that, at 
low densities, the derivative of the energy with respect 
to density is determined by the spatial average of the 
two-body interaction, and this changes from positive to 
negative as a increases through ac2 ■ The contribution to 
the kinetic energy varies as for small and is there- 
fore unimportant for weak coupling. For fermions, the 
energy per particle is always positive for small because 
the Fermi energy dominates the total energy. 

Correlations — We now explore the nature of the 
ground state in situations where particles are concen- 
trated near 9 = ±7r/2. Classically, for a repulsive in- 
teraction one would expect essentially equal numbers of 
particles to be present on each side of the ring, while for 
an attractive interaction, the particles would be concen- 
trated on one side of the ring. A convenient diagnostic 
for probing the nature of the state is the pair distribution 

function p2 (6', 6*') = T,^^j{H0-^i)H0' -^j))- For weakly 
interacting bosons and < a < ac, when the interac- 
tion is positive and homogeneous, the system is well de- 
scribed as a Bose-Einstein condensate, with all particles 
in a single-particle state that has density maxima at the 
points 9 = ±7r/2. As shown in Fig. [BJ the pair distribu- 
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Figure 6: Pair distribution function ^2(^1,^2) for A*' = 4 
bosons at O! = 0.207r and nrd = 0.026 (left), nrd = 0.803 
(center) and nrd ~ 0.808 (rigiit). Black corresponds to p2 = 
and white to the highest values. 

tion function has maxima in the vicinity of 6, 6' — ±7r/2, 
and they are all of equal strength. For stronger coupling, 
tunneling between 9 — ±7r/2 is suppressed, and atoms 
become localized in the vicinity oi 9 = ±7r/2. In this 
regime, it is more appropriate to think in terms of Mott- 
insulator-type states |iV+, A^_) with defiirite numbers of 
particles N± localized near 9 — ±tt/2. The ground state 
for a repulsive interaction would be expected to be dom- 
inated by I-/V/2, N/2), where for simplicity we have taken 
N to be eveir. In this case, oire expects the peaks arouird 
6* « 6*' sa ±7r/2 to have a strength 1-2/iV (1/2 for iV = 4) 
times that of the peaks a.t 9 ^ —9' . This describes well 
the data for nrd = 0.803. 

Particles in a bound state would be expected to local- 
ize oir oire side of the ring, iir the states |A^, 0) or |0,iV), 
in order to maximize the attractive contribution to the 
energy. Because of the symmetry of the system under 
replacement of 9 by —9, the ground state will be a lin- 
ear combiiration of these states with equal probabilities. 
For such a state, the pair distribution function is qualita- 
tively different from that of unbound states iir that it has 
only two peaks, at 9 ^ 9' ±7r/2, as is illustrated in Fig. 
inifor nrd ~ 0.808, a value only slightly different from the 
previous one. For very strong coupling nrd 3> 1, the crit- 
ical angle for the formation of bound states approaches 
ac as is seen in Fig. For weaker coupling, the critical 
angle for bound state formation iircreases because of the 
limited range of CM angles around 8 = ±7r/2 for which 
the interaction is attractive. 

Finally, we discuss some experimental aspects. Quasi- 
ID ring traps with R/a± ~ 10 — 10^ and a_L ~ 1 — lO^m 
may be produced using techniques such as time-averaged 
optical tweezers [lo| or radiofrequency-dressed magnetic 
traps (ll| . The weak-coupling regime nrd ^ 1 may 
be explored with magnetic atoms (for Cr, ^ Inm), 
but to observe strong correlations (7 = TiRrd/Na\ ^ 
1) it is preferable to work with only a small number 
{N ^ 10 — 100) of molecules and to increase R/a±. The 
strongly coupled regime nrd 3> 1 may be realized with 
polar molecules (r^ ^ 0.1/im for KRb) and large parti- 
cle numbers N ^ 10^i?//im (at a = 0). Measurement 
of clock shifts similar to those done to determine the 
atom numbers on sites in the superfluid-insulator tran- 



sition [18| would be a useful way to distinguish between 
the state \N/2, N/2) and Schrodinger-cat-like states that 
are superpositions of l-A'^, 0) and |0, -/V). 

In conclusion, we have shown that dipolar particles 
in ring traps may be used to realize a variety of dif- 
ferent states due to the interplay between the ring ge- 
ometry, dipolar anisotropy, and confinement. This in- 
cludes both inhomogeneously localized states as well as 
unconventional effective short-range physics, which may 
be explored by tuning the numbers of dipolar atoms or 
molecules, their orientation, and the trapping parame- 
ters. 
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